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Abstract 

. We describe a simple algorithm that calculates the distributions of electrons 

and positrons under multiphoton beamstrahlung at a future linear collider. The 
fSj ' evolution equation as given by Chen is solved by a Monte Carlo algorithm. Explicit 

• multiple beamstrahlung photons are generated. We present first results from an 

implementation of beamstrahlung into the Monte Carlo event generator WOPPER, 
that calculates the QED radiative corrections to the process e~^e~ — > 4 fermions 



p • through resonating W pairs. 

!• 1 Introduction 



At future linear e~^e~ colliders with center of mass energies of 0.5 TeV and above, beam- 
strahlung |1|, 0, the synchrotron radiation from the colliding e~^e~ beams, will become of 
significant physical interest. On one hand, beamstrahlung may carry away a substantial 
fraction of primary beam energy and lead to a degradation of the effective center of mass 
energy for e~^e~ collisions, on the other hand, the lower energy e^e~ and 7's contribute 
to background processes. 

In general, the full treatment of beam-beam interactions and beamstrahlung is a com- 
plicated many-body problem, and dedicated codes for the calculation of beam-beam effects 
exist [0] or are being developed |Q. However, from the high energy physics perspective 
the availability of a simplified analytical or numerical treatment is desired for the study of 
interesting processes, (e.g., e~^e~ — >■ W^W~), at the level of Monte Carlo event generators. 

* Contribution to the Proceedings of the Workshop on Physics with e^e^ Linear Colliders, Annecy, 
Gran Sasso, Hamburg, February 4 to September 1, 1995. 
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As long as the average number of beamstrahlung photons per beam particle is much 
less than unity, the radiation spectrum is given by the Sokolov-Ternov spectrum 0. In 
the case of the proposed linear colliders, the average number of radiated photons per beam 
particle at design luminosity is typically 0{1). Therefore, the effects of multiple radiation 
of beamstrahlung photons have to be taken into account. 

Multiphoton beamstrahlung has already been discussed in ^ (and references quoted 
therein). In this note we will follow the approach by Chen et al., however, without using 
any approximation to the radiation spectrum. We shall derive a recursive solution of the 
evolution equation, which is well suited for a numerical implementation. 

The Monte Carlo algorithm for beamstrahlung presented below is considered comple- 
mentary to a full treatment of all beam-beam effects as performed in 0. The main purpose 
of this Monte Carlo approach is to provide a simple package that can be incorporated into 
or interfaced with Monte Carlo event generators for physics at future e+e~ colliders. In 
section 2 we will review the electron energy spectrum under multiphoton beamstrahlung. 
Section 3 outlines our Monte Carlo approach, while in section 4 we present preliminary 
results. Further details of the calculation will be presented elsewhere [||. 



2 Electron Energy Spectrum 



In this section we recapitulate the derivation of the electron energy spectrum due to 
multiphoton beamstrahlung in the spirit of Chen et al. 0. We use h = c = 1. 

Let ip{x, t) be the distribution function of the electron with energy fraction x = E/Eq 
at time t, normalized such that 

dxilj{x,t) = 1 . (1) 

Jo 

We shall assume that the emission of the photons takes place in an infinitesimally short 
time interval. Therefore the interference between successive emissions may be neglected, 
and one can describe the time evolution of the electron distribution by the rate equation 

^ = -iy{x)tlj{x, t) + r dx' F{x, x')tlj{x', t) , (2) 

ut Jx 

where the first and second term on the r.h.s. correspond to the sink and source for the 
evolution of ip, respectively. F is the spectral function of radiation, and 

u(x) = r dx" F(x",x) (3) 
Jo 



2 



represents the average number of photons radiated per unit time by an electron with en- 
ergy fraction x. Electrons and positrons from pair production in photon-photon processes 
are not taken into account. 

The spectral function of radiation is conventionally characterized by the beamstrahlung 
parameter T, which is defined as 

nie 

where B is the effective field strength in the beam, and Be = m^/e is the Schwinger 
critical field. The e~^e~ beams generally have Gaussian charge distributions, and thus the 
local field strength varies inside the beam volume. However, it has been shown through 
integration over the impact parameter and the longitudinal variations, that the overall 
beamstrahlung effect can be simply described as if all particles experience a uniform mean 
field Bmean duriug an effective collision time r = leff/'^ = V^o'z- In what follows, we will 
assume the obtained mean beamstrahlung parameter T^ean and effective beam length 
leff but drop the subscripts. 

The spectral function of synchrotron radiation was derived by Sokolov and Ternov |^ 
and reads 

where ^ = (3T/2)x', r] = (2/3T)(l/x — and K^, denotes the modified Bessel function 
of order u. 

We define 

iy{x) = / dx' F{x',x) = Uci ■ Uo{xT) , (6) 
Jo 

where Ud is the number of photons per unit time in the classical limit, 

5 ami 

which is independent of the particle energy for a given field strength due to (^). 
The function Uo{y) behaves as 



Uo{y) 



1 y^O 
const. X y~^^^ y ^ oo 



exhibiting the asymptotic approach to the classical limit and the suppression of radiation 
in the deep quantum regime, respectively. 
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Finally, the differential e~^e luminosity can be expressed as the convolution of the 
effective electron energy distributions of the colliding beams 0, 

dC(s) f 

= ^0 J dXidX2 5{s - XiX2So)'4){xi)%l){x2) , (9) 

where Cq is the nominal luminosity of the collider, including the enhancement factor due 
to beam disruption, and the effective energy distribution il){x) is obtained by averaging 
over the longitudinal position within the beam, 

2 /-Vs 

iIj{x) = - ijjix^t). (10) 



/ Jo 



3 Multiphoton Beamstrahlung by Monte Carlo 

The evolution equation (^ may be solved by a Monte Carlo algorithm. Here we will 
outline only the basic idea, deferring the full presentation and discussion to a forthcoming 
pubhcation p. 

To determine the electron distribution ilj{x,t), let us decompose it into a suitable set 
of time independent functions 0ri(x) using the ansatz 

^{x,t) = J2Mx)--^e-''\ (11) 

n=0 "'• 

where i/ is a parameter that will be determined below. The normalization condition (|T]) 
translates into 

1^ dx (t)n{x) = I . (12) 

Jo 

Plugging the ansatz ( pT]) into the evolution equation (H) we obtain a recursion relation 
for the (/)„, 



0„(x) = / dx' 
Jo 



v{x)\ ^, ,, F(x,x') 
1 ^ J S{x - x') + ^ J ^ 



'M-l 



'x') for n > 1 , (13) 



with initial condition 

(/)o(x) = ^(x,0). (14) 

If we neglect the intrinsic (Gaussian) energy spread of the accelerator, the initial condition 
of the electron distribution is simply given by 

V'(x,0) = 5(1 -x) , (15) 
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and thus 



<Po{x) = 6{1 — x) 



(16) 



The functions (j)n{x),n > 0, can be determined, e.g., by numerical integration. Here 
we suggest a different way to exploit the recursion relations. First note that the kernel 
given by the square brackets in ([13|) is positive for x' > x. Furthermore, since z/(0) > 
z/(x) > the kernel is positive if we choose u = z/(0) = u^. The functions 0„ will 

then be positive, and together with the normalization condition (|T^ ) they will allow for a 
probabilistic interpretation. 

Using our ansatz (p!T|), the effective electron energy distribution (|T0|) reads 

-f{l + n;Nciy 



w. 

where 7(1 + n; Nd) is the incomplete gamma function, and A*",,; = t'ci ■ is the average 
number of photons per beam particle radiated during the entire collision, in the classical 
limit. Due to the positivity of the term in square brackets in ([l7D and the functions 0„, 
the electron energy distribution can be determined straightforwardly by a multichannel 
Monte Carlo method M. 



4 Results 

Of the abundance of available parameter sets for future linear e~^e~ colliders, we have 
taken a hopefully representative selection from |^ that comprises many of the proposed 
linear collider designs. Table ^ lists some of the parameters that are relevant for our 
numerical results. Tmean was determined similarly to 0. 

We have calculated the differential luminosity of the e^e^ beams under beamstrahlung, 
eq. (^, using our Monte Carlo algorithm based on (|l^). The result of the simulation for the 
normalized luminosity is shown in figure [l|. For the sake of brevity we have selected only a 
subset of four designs which essentially cover the interesting parameter range: the TESLA 
design with rather large bunch length but small effective beamstrahlung parameter and 
a correspondingly rather soft beamstrahlung spectrum, the JLC-S design with shorter 
bunches but large T and large mean energy loss due to a hard beamstrahlung spectrum, 
NLC with short bunches and moderate T, having a small number of radiated photons A*",,;, 
and VLEPP with a large multiplicity of beamstrahlung photons. We have also calculated 
the relative luminosity for collisions with s > 0.99 Sq, the result being displayed in the 
last line of table |l[ 
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Parameter 


SBLC 


TESLA 


NLC 


JLC-S 


JLC-X 


CLIC 


VLEPP 


Luminosity [10^'^cm~^s^^] 


5.3 


6.1 


6.9 


4.6 


5.1 


? 


9.7 


Particles/bunch [10^°] 


1.1 


3.63 


0.65 


1.44 


0.63 


0.8 


20 


ax [nm] 


335 


845 


320 


260 


260 


250 


2000 


ay [nm] 


15.1 


18.9 


3.2 


3.04 


3.04 


7.5 


3.9 




300 


700 


100 


120 


90 


200 


750 


T 

mean 


0.049 


0.029 


0.090 


0.213 


0.120 


0.072 


0.064 


Na = Vci ■ 1/2 


1.42 


1.93 


0.87 


2.47 


1.04 


1.4 


4.6 


C{s > 0.99 So) /Co 


0.58 


0.55 


0.68 


0.32 


0.62 


0.56 


0.18 



Table 1: Selected parameters of some proposed linear collider designs. 



In high energy e~^e~ processes, besides beamstrahlung one also has to take into ac- 
count the radiative corrections to the process under consideration. The single most 
important universal contribution of the radiative corrections are the leading logarithms 
a" log"(s/mg), which are due to QED initial state radiation (ISR). These large logarithms 
can be resummed to all orders in the structure function approach, with the electron struc- 



ture function 10 



D{z; /i^) = §(1 - z)^-' (l + - + z) + 0{(3') , (18) 

where 

;3 = H^flog4-l) . (19) 

/i^ is the factorization scale of order s, and z is the electron energy fraction after initial 
state radiation. 

As an application, we have interfaced our Monte Carlo routines for beamstrahlung 
with the Monte Carlo event generator WOPPER [11|. WOPPER simulates the process e'^e~ —>■ 



(W^W~) —>■ 4/ including QED radiative corrections in the abovementioned leading loga- 
rithmic approximation with explicit photons. In figure ^ we show the normalized invariant 
mass distribution M^^y of the final state (assuming a 100% reconstruction efficiency for 
simplicity) , 

— = / dxidx2dzidz2 i'{xi)i){x2) D{zi\ ^^)D{z2\ ^i^) ■ —ao{xiZiX2Z2s) , (20) 
ao dM^^ J ao 

using the same collider designs as in figure |l|. The thin continuous line represents the 
corresponding distribution from ISR alone. For events with a significantly reduced in- 
variant mass the corrections due to beamstrahlung are numerically much larger than due 
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e'^e Luminosity at 500 GeV 

JLC-S 




Figure 1: Differential e'^e luminosity. 



to pure ISR, and the shape of the distribution is essentially dominated by the beam- 
strahlung spectrum. This result clearly shows the importance of beamstrahlung and the 
need to include beamstrahlung into Monte Carlo event generators for physics at future 
e'^e~ colliders. 
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